Holonomic quantum computation in decoherence-free subspaces 
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We show how to realize, by means of non-abelian quantum holonomies, a set of universal quantum 
gates acting on decoherence-free subspaces and subsystems. In this manner we bring together the 
quantum coherence stabilization virtues of decoherence-free subspaces and the fault-tolerance of 
all-geometric holonomic control. We discuss the implementation of this scheme in the context of 
quantum information processing using trapped ions and quantum dots. 



PACS numbers: 03.67.Lx,74.20.Fg 

Introduction. — The implementation of quantum infor- 
mation processing (QIP) poses an unprecedented chal- 
lenge to our capabilities of controlling the dynamics of 
quantum systems. The challenge is twofold and some- 
what contradictory. On the one hand one must (i) main- 
tain as much as possible the isolation of the computing 
degrees of freedom from the environment, in order to pre- 
serve their "quantumness" ; on the other hand (ii) their 
dynamical evolution must be enacted with extreme pre- 
cision in order to avoid errors whose propagation would 
quickly spoil the whole quantum computational process. 
To cope with the decoherence problem (i), active strate- 
gies such as quantum error correcting codes |]J, as well 
as passive ones such as error avoiding codes |2(, have 
been contrived. The latter are based on the symmetry 
structure of the system-environment interaction, which 
under certain circumstances allows for the existence of 
decoherence-free subspaces (DFS), i.e., subspaces of the 
system Hilbert state-space over which the dynamics is 
still unitary. DFSs have been experimentally demon- 
strated in a host of physical systems (e.g., 0,0]) The 
DFS idea of symmetry-aided protection has been gener- 
alized to noiseless subsystems [f|, experimentally tested 
in Ref. 0. 

Holonomic quantum computation (HQC) [?j is an all- 
geometric strategy wherein QIP is realized by means of 
adiabatic non-abelian quantum holonomies Quan- 
tum information is encoded in a degenerate eigenstate 
of a Hamiltonian depending on a set of controllable pa- 
rameters e.g., external laser fields. When the latter are 
adiabatically driven along a suitable closed path, the ini- 
tial quantum state is transformed by a non-trivial unitary 
transformation (holonomy) that is geometrical in nature. 
Following the ion traps HQC-implementation proposal 
, several other schemes, based on a variety of physical 
setups, were proposed |lCj . One expects the geometri- 
cal nature of quantum holonomies to endow HQC with 
inherent robustness against certain errors. This alleged 
fault-tolerance has only recently been subjected to seri- 



ous scrutiny the resulting, still developing picture, is 
that while stability against decoherence must be further 
assessed (indeed the adiabatic theorem was only recently 
generalized to open quantum systems [12|), HQC seems 
to exhibit a strong robustness against stochastic errors 
in the control process generating the required adiabatic 
loops 13]. From this point of view HQC seems to be 
promising with respect to the general challenge (ii). 

In this work we describe a QIP scheme which combines 
DFSs and HQC Hi \u&. More specifically, we show how 
to perform universal quantum computation within a two- 
qubit DFS for collective dephasing by using non-abelian 
holonomies only. The discussion is then extended to con- 
sider general collective decoherence as well. The appeal 
of such a strategy should, in view of the above, be evi- 
dent: try to bring together the best of two worlds, namely 
the resilience of the DFS approach against environment- 
induced decoherence and the operational robustness of 
HQC. Moreover, we formulate our results using rather- 
generic Hamiltonians, so that the scheme proposed in 
this work appears to be a suitable candidate for experi- 
mental demonstration in a variety of systems, including 
trapped ions and quantum dots. 

Dark-states in a Decoherence- Free subspace. — Let us 
start by considering a four-qubit system with state- 
space Hi = (C 2 )® 4 . We denote by X U Y X and Zi 
the three Pauli matrices acting on the /th qubit; for 
any pair 1 / ra of qubits we define the operators 
R?m := 5 {XiX m + YiY m ), R v lm :— % (XiY m — YiX m ), 
Rf m :— I (Z m — Zi). These operators have a non-trivial 
action over the subspace of Hi m — C 2 ® C 2 ^ spanned by 
|0/l m ) and |l;O m ) (where |0)/|1) arc the +/— eigenstates 
of Z, respectively): they provide a faithful representa- 
tion of the su(2) algebra over this subspace, where they 
act as the Pauli matrices, while they vanish on the or- 
thogonal complement spanned by |0;O m ) and |l;l m ). Let 
Z := £i=i Z h then [Rf m ,Z] = (V/, m). It follows that 
every eigenspace of Z is invariant under the action of the 
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R^'s. In particular this holds true for the subspace 

C := span{|1000), |0100), |0010), |0001)}. (1) 

C is a DFS against collective dephasing 0], i.e., states in 
C are immune from decoherence induced by system-bath 
interactions of the form Z £g) £?, where B is an arbitrary 
bath operator. Collective dephasing is known to be a 
major source of decoherence in ion- trap based QIP Q ■ 
We assume that the system dynamics is generated by 
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where m are qubit indices and the J; m 's are controllable 
coupling constants. These are the parameters that will be 
driven along controlled adiabatic loops to enact quantum 
gates via non-abelian holonomies. When Jf = 0, H 
is the XY Hamiltonian found in a variety of quantum 
computing proposals, e.g., the quantum Hall proposal 
[lif . quantum dots ^3 an d atoms in cavities |lSj . It also 
describes trapped ions subject to the S0rensen-M0lmer 
scheme 0. The case Jf m = is related to the XY 
model via a unitary transformation. 

By way of introduction to our HQC-DFS scheme, note 
that Hamiltonian J2J) has a hidden multi-level structure 
with interesting properties (isomorphic to the one ex- 
ploited in Ref. 9]). Indeed, the Hamiltonian Hi„ m — 

Y,j=m n J jl( R ji cos fji + R % s ' ml P3i), in the basis {|e) := 
|100) 

Irani | ffl 

) :— |010)z m „, I02) |001); mn }, with I ^ 
m =/= n, takes the form 



H lmn = Ji m e iv '™\e)( gi \ + Ji n e i ^\e){g 2 \ +h.. 
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This is a so-called Lambda scheme, with |e) at the top 
and |<7i), \g 2 ) at the bottom. It is well known that for 
every value of the Jji's and fjis, Hamiltonian Q has one 
dark state \D(J lm ,J ln )) oc Ji m e i,pim \g 2 ) - 4^'° \9i), 
i.e., a state satisfying H lmn \D(Ji m , J ln )) = 0. 

The key idea is now as follows: moving in the control 
parameter space to the point Ji m /Jin = 0, one has that 
the dark state is given by \gi); then, if one adiabatically 
changes the parameters along a closed loop, at the end 
the state \gi) will pick up a non-trivial geometric phase. 
This simple fact, suitably generalized, is the basic ingre- 
dient of our HQC scheme. We further supplement the 
dark state |<7i) by another state, that is also annihilated 
by the system Hamiltonian (and thus does not acquire a 
dynamical phase either); these states together will form 
our qubit. In other words, the crucial observation is that 
one can embed within the DFS C [Eq. 0J] a manifold of 
dark states that, for specific values of the coupling con- 
stants in Eq. J3J), coincide with logical encoded qubits. 
Holonomic manipulations are then used to generate a 
universal set of quantum logic gates. 

Let us stress that even if during state manipulation the 
system described by Eq. leaks out of the logical en- 
coding subspace (the states \0)l and defined below), 



e.g., due to the breakdown of adiabaticity [EJ, this re- 
sults just in the kind of errors that HQC is robust against 
Moreover, during such leakage the system never 
abandons the DFS C. Thus protection against collective 
dephasing is maintained throughout the whole gating pe- 
riod, which in turns allows one to stretch the gating time 
in such a way as to fulfill the adiabatic constraint for 
a longer period of time than would be possible without 
the DFS. These remarks at least partially counter a stan- 
dard objection to HQC, that the use of slow gates gives 
decoherence more time to exert its detrimental effects. 

One-qubit holonomic gates. — We now show how to 
realize the single-qubit phase gate exp(iipZL), where 
Z L \a) = {-l) a \a) L (note that Z L = Rf 2 ). We encode a 
logical qubit in C by \1) L := |0010) and |0) L := |0001), 
while the remaining pair |oi) := 1 1000) and \a 2 ) := |0100) 
plays the role of ancillae. Let us set all the Ji m to zero, 
except J24 and J34, such that the Hamiltonian J5J) re- 
duces to 

H z = J 2 a(R 24 cos Lp 24 + R y 2i sin (p 24 ) + J Z4 R X M . (4) 

We can also write Hz = -/34|ai)(a2|+J r 24e lv |ai)(l|L+h.c. 
This is a Lambda configuration with \a\) at the top and 
|o2))|1)l at the bottom. Therefore, as in our discus- 
sion above, Hz has a zero-eigenvalue eigenstate (dark 
state) given by = cos6*|1)l — sin^e^le^) where 

9 — tan _1 (J 2 4/J34) and <p = ip 24 . The state |0)l is also 
a zero-eigenvalue eigenstate that does not depend on the 
parameters 9 and tp. By adiabatically changing cp in such 
a way as to have a loop starting from 9 = ip = 0, the state 
acquires a Berry phase which is proportional to the 
solid angle Clz swept out by the vector (0, 9) [2(|. There- 
fore after this adiabatic loop one has that: |0)l — > |0)i 
and — > e~ lflz ^ 2 |1) l, which is clearly equivalent to 
the operation exp(iri^Zi/2). 

In order to obtain a universal set of gates we need to 
generate at least two non-commuting single-qubit gates. 
Therefore we next show how to implement exp^ipX^) = 
cxp(itpRf 2 )- One way is to once again establish an iso- 
morphism between the system governed by Eq. J5J), re- 
stricted to C, and the HQC model of Ref. [9j. Here we 
provide an independent derivation. We turn on the cou- 
plings in such a way as to obtain the Hamiltonian 



H X = </34i?34 + J 2 4{C0S ifi 



H 24 



V2 



^14 + sin y ^4 ^14 y 



V2 



( 5 ) 

Let \±) L := (\1) L ± \0) L )/y/2. Then one can readily 
check that under the action of (R 24 — i?f 4 )/v / 2 the state 
i-> 0, \-) L 1-* |oi), and |ai) 1-* \-)l, etc. There- 
fore Hx = J34|ai)(a 2 | + J 24 e iip \ai)(— \l + h.c. This is a 
Lambda configuration with \a%) at the top and |a 2 ), |— ) j, 
at the bottom, so that Hx supports a dark state l^) = 
cos0|-)i - sinfle^lc^) where 9 = tan -1 (J 24 / ^34) . The 
similarity between Hx and Hz is evident. Then, by exe- 
cuting an adiabatic loop in the parameter space in anal- 
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ogy to the Hz case, one obtains the geometric evolution 
-> \+)l, Hl -» e- l ^/ 2 |-) L , where now ft x is 
the solid angle swept out by the vector (9,ip). Switch- 
ing back to the computational basis this transformation 
amounts to the map |0)l — * cos^plO)^ + isin^ L |l)i, 
and —* cos^-|l)x + ism^2L\0) L , which is equiva- 
lent to exp(if2x^i2/4)- 

Two-Qubit holonomic gates. — A crucial, and typically 
rather demanding part of any QIP implementation pro- 
posal, is the realization of an entangling two-qubit gate. 
We next consider this problem and demonstrate how to 
solve it following a strategy relying on the same abstract 
holonomic structure as that discussed above for one-qubit 
gates. The total state-space is now Ttf 2 , while the two- 
qubit code is spanned by \a) l ® \0)l, (a, = 0, 1). The 
states |ai)® 2 (i = 1,2) will function as ancillae. Impor- 
tantly, once again all the relevant states belong to a DFS 
against collective dephasing. 

Let us suppose that one can engineer the following con- 
trollable four-qubit interaction 

Hi = J24 : 68(^24 cos f + -^24 sin f)(^6S cos <P + R 6S sin <P) 
+ ^34,78-^34^78' (6) 

which should be recognizable as a straightforward exten- 
sion of the one-qubit Hamiltonian (J2J. Below we dis- 
cuss the implementation of such a Hamiltonian. To ex- 
plicitly exhibit the dark state structure, we write this 
Hamiltonian in the form H4 = J 3 4,78-R34(|ai)(a 2 |)® 2 + 
^24,68e IV (|ai)(lL|)' 812 + h.c, from which it is easily seen 
that there is one dark state, given by cos#|l)® — 
sin6»e iV |a 2 )® 2 , where 6 = tan" 1 (J 24 , 6 8 / 'Jm,7s) • After the 
same kind of adiabatic cyclic evolution as described in 
the single-qubit gate case, one obtains in C® 2 the con- 
trolled phase-shift gate CP = diag(l, 1, 1, e lfip / 2 ), where, 
as usual, flp is the solid angle swept out by the vector 

Extensions and generalizations. — The scheme de- 
scribed so far for the case of collective dephasing is 
straightforwardly scalable to an arbitrary number TV 
of encoded DF qubits. The total space is now given 
by C 8,v = (C 4 )® N , with C given in Eq. (JTJ. This, 
of course, is still an eigenspace of the collective spin 
^-component, i.e., Z = Ya=i ^i- Following the pro- 
cedure established above, the controllable Hamiltonian 
used to generate a controlled phase-shift between the 
j-th and the j-th encoded qubits, has the same struc- 
ture as H4 [Eq. ©], where, with obvious notation, 

>^24;68 - * ^4(i-l)+2,4(i-l)+4:4(j-l)+2,4(j-l)+4 and J34;78 ~* 

^4(i-i)+3,4(- i -i)+4;4(j-i)+3,4(j-i)+4, and similarly for the 
RkliR-kl operators. 

Next, by making use of noiseless subsystems (NS) 
13 , we show that our combined HQC-DFS strategy can 
also be applied against general collective decoherence . 
Our arguments will be existential in nature, with con- 
structive details to be discussed elsewhere. Let TCj = 



Q2.7+1 denotg the total spin- J irreducible representation 
of su(2). The state-space of five qubits, i.e., Hf, 2 , decom- 
poses, with respect to the collective sw(2)-representation 
(Clebsch-Gordan decomposition) as follows 

nf f l^ c 4 ®w 3 /2©c 5 ®Hi/2©c®h 5/ 2. (7) 

Each of the C nj factors represents the multiplicity of the 
total spin- J irreducible representation and corresponds to 
a NS against collective general decoherence [f| . Consider 
the first term in 10: the multiplicity factor C 4 for the 
J = 3/2 representation provides a four-dimensional NS. 
It might then encode two noiseless qubits, but, since we 
wish to perform QIP with holonomies, we will instead use 
this C 4 space as a code for just one noiseless qubit \ch)l, 
(a = 0, 1) and two ancillary states |5,-), (i = 1,2). Sup- 
pose now that one is able to enact the controllable Hamil- 
tonian H m = J'|ai)(a 2 |+J^|ai)(0| L +J( / e^|ai){l| i + 
h.c, which, when Jq — 0, admits a dark state given by 
|*3) =cos0|l)i -sin0e^|o2>, where 9 = tan~ 4 ( J'{/ J'). 
By resorting to the same considerations as previously de- 
veloped for the collective dephasing case, it should be 
clear that this allows us to enact a phase gate Zl be- 
tween |0)l and By choosing Jq,J" such that H^s 
becomes J'|oi)(a2| + J'le l(p \d>i)(—\L + h.c. we can enact 
the Xl gate, so that universal single-qubit control by 
holonomies can be achieved in this case as well. 

To realize the required multi-level controllable Hamil- 
tonian, we observe that the C 4 space under considera- 
tion is a four-dimensional irreducible representation of 
the permutation group S$ [acting over the whole space 
as a <Sij —i \5 u) = ®] = i\a a {j)) , wherecr € S5]. Follow- 
ing Ref. [21|, universal control over this irreducible rep- 
resentation space amounts to the ability to switch on 
and off a pair of generic Hamiltonians in the group 
algebra of S5. An important example is provided by 
Heisenberg exchange Hamiltonians, i.e., X)/<m Jim^r^m 
[where S; = (JQ, Yi, Zi)/2], which are naturally avail- 
able interactions in several spin-based proposals for QIP 
[22I E^ . The construction of two-encoded-qubit holo- 
nomic gates, and the generalization to arbitrary numbers 
of such qubits, again follow the same pattern as in the 
collective dephasing case. 

Implementation. — We note that all required Hamil- 
tonians, Hz, Hx and H^, have a similar form. These 
Hamiltonians all involve control over both 9 and ip. How- 
ever, we are free to choose any loop C{9, ip) in the {9, <p)- 
parameter space, and we can choose a loop which tog- 
gles between 9 and p. For example, the loop C(9, ip) : 
(0,0Mf,0Mf,<A)M0,<A)) = (0,0) has this prop- 
erty, where <po is the solid angle f2 swept out by the vector 
(9, ip) for this specific loop. 

Let us briefly address the feasibility of the control 
of Hz, Hx and H4 in the context of actual quantum 
computing proposals. In proposals based on electron 
spins in quantum dots |22|, the Hz and Hx Hamilto- 
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mans are available when one takes into account the spin- 
orbit interaction, for then the effective spin-spin interac- 
tion becomes ffy(t) = JiAt)\B(t)(X 4 Yj - YiXj) + (1 + 
^(t^iXiXj+YiY^+ZiZj] OUg. This Hamiltonian has 
enough degrees of freedom to implement Hz and Hx , via 
separate control of J and the dimensionless anisotropy 
parameters (3 and 7, assuming the latter can be made 
^ 1, e.g., when the inter-dot distance is large [24]. The 
parameters (3 and 7 can be controlled via the confining 
potential or viapulse shaping, as has been discussed in 
detail in Ref. |2g, and ip — tan _1 /9/(l + 7). The pa- 
rameter 9 is controllable via, e.g., J24 and J34. Imple- 
mentation of Hi is undoubtedly more challenging, since 
four-body interactions are involved. Recent results in- 
dicate that such terms arise by simultaneously coupling 
four quantum dots pif: their relative strength can be 
controlled, e.g., by adjusting the confining potentials [2^ . 

In proposals based on trapped ions, the implemen- 
tation of Hz and Hx is directly possible using the 
S0rensen-M0lmer (SM) scheme [13 . The SM gate be- 
tween two ions implements Hz and Hx with control over 
the various terms achieved via the phases of two lasers 
[28|. It is also possible to implement H4 u sing the SM 
scheme, via control over two pairs of ions [28j. Given 
that a geometric two ion-qubit phase gate has already 
been demonstrated |29j | , trapped ions seem to be partic- 
ularly favorable for the implementation of our proposed 
HQC-DFS scheme. 

Conclusions. — We have combined the DFS and HQC 
techniques, and showed how to realize universal quan- 
tum computation over a scalable DFS against collective 
dephasing by using adiabatic holonomies only. We dis- 
cussed an extension to the general collective decoherence 
case, arguing that controllability of exchange Hamiltoni- 
ans would suffice. Remarkably, the whole computational 
process is carried out within the DFS. The DFS em- 
bedding, along with the all-geometrical nature of HQC, 
promises to give to this scheme a twofold resilience, 
against decoherence and stochastic control errors. The 
proposed universal quantum gates are carried out via 
adiabatic manipulation of Hamiltonians that are control- 
lable in several proposals for QIP implementations. We 
are therefore hopeful that the theoretical ideas presented 
here may stimulate corresponding experimental activity. 
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